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Abstract 



We consider the nonleptonic B decays B -> £>W.D S (2317) and B ->■ L>MD S (2460), involving the 
newly discovered Z? s (2317) and the Z3 S (2460) states. We find that experiments indicate disagree- 
ment with model calculations of their properties and/or breakdown of the factorization assumption 
for these decays . We point out that decays involving B s mesons where the D s resonances can be 
produced via the weak decay of the b quark can provide further information about the nature of 
these newly discovered states. We also propose a model to calculate the two body nonleptonic de- 
cays B -> £>WL> S (2317)(L> S (2460)), if the D s (2317) and £> s (2460) are interpreted as DK and D*K 
molecules. 



1 Introduction 



There has been recent observations of an unexpectedly light narrow resonance in D+tt° with a mass of 
2317MeV/cP by the BaBar collaboration [1], together with another second narrow resonance in DgTr ^ 
with a mass 2460MeV/c 2 [2]. 

The smaller than expected masses and narrow widths of these states have led, among other explanations 
[3] , to a multi-quark anti-quark or a DK molecule interpretation of these states [4] , or to an interpretation 
as p-wave states where the light degrees of freedom are in an angular momentum state j q = \ [5], or 
even some combination of these [6]. There are also conflicting lattice interpretations of these states [7]. 
The mass difference between the D s (2317) and the well established lightest charm-strange meson, D s , is 
AM = 350MeV/c 2 . This is less than the kaon mass, thus kinematically forbidding the decay A, (2317) — > 
D u ,d + K. The possible resonance at 2460MeV/c 2 also has such a mass difference when taken with the 
lighter D* state. The interpretation of these states as bound D^K molecules just below the D^K 
threshold is particularly interesting in the light of the recent discovery of a narrow resonance in the decay 
J/tp — > 7pp [8] which has been interpreted as a zero baryon number, "deuteron-like singlet 1 <S'o" bound 
state of p and p [9] . 

In the heavy quark theory, the ground state heavy meson involving a heavy and a light quark has the 
light degrees of freedom in a spin-parity state jq = \ , corresponding to the usual pseudoscalar-vector 
meson doublet with J p = (CP, l - ). The first excited states involves a p-wave excitation, in which the light 
degrees of freedom have j p = \ + or | + . This leads to two heavy doublets, the first giving J p = (0 + , 1 + ) 
and the latter a heavy doublet with J p = (l + ,2 + ). Heavy quark symmetry rules out any pseudoscalar 
coupling of this doublet to the ground state at lowest order in the chiral expansion [10] and so these states 
are expected to be narrow. Recent Belle analysis of B~ — > D^*} n + ir~ decays [11] indicate the presence 
of the 1+ state in this multiplet at a mass of M D o = (2421.4 ± 2.0 ± 0.4 ± 0.8)MeV/c 2 with a width of 
F D o = (23.7 ± 2.7 ± 0.2 ± 4.0) MeV . The other state in the doublet (2+) is also found with a mass of 
M D o = (2461.6 ± 2.1 ± 0.5 ± 3.3)MeV/c 2 with a width of T D o = (45.6 ± 4.4 ± 6.5 ± 1.6) MeV . In the 
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D s system the counterpart states to these are naively expected to be a 100 MeV heavier because of the 
strange quark mass and so these states can probably be identified with D s i(2536) and D s j(2573) [12]. This 
is in line with the experimental observations that in the ground state the D s mesons are about a 100 MeV 
heavier than their nonstange counterparts. 

The other excited doublet has J p = (0 + ,l + ). These states are expected to decay rapidly through 
s-wave pion emission in the & system and by kaon emission in the D s system and have large widths 
[13]. Observation of the 1 + state in the D system was reported by CLEO [14] some time ago. The 
recent Belle analysis of B~ — > D ( - + *\ + ir" decays [11] also find evidence for the states in this doublet at 
M D »o(0+) = (2308 ± 17 ± 15 ± 20)MeV/c 2 with a width of F D o = (276 ± 21 ± 18 ± 60) MeV . The other 
state in the doublet is also found with a mass of M D *o(l + ) = (2427 ± 26 ± 20 ± 15) MeV/c 2 with a width 
of F D o = (384i^° 7 ± 24 ± 70) MeV . Note that these states are broad as expected from theory. Naively 
then, we should expect the D s counterparts of these states at Mn s (0 + ) rs 2408 and Md s (1 + ) k, 2527. 
These numbers are consistent with quark model estimates [15] and we expect these states to be broad. 

The recently observed D s resonances have masses below these expectations and are very narrow, decaying 

(*) 

through isospin violating transitions to Dl 'ir final states. This has generated speculations that these states 
may not be p-wave excited states but rather something exotic like D^K molecules. 

While the spectroscopy of these newly discovered states can provide clues to their structure, decays 
involving these states can yield further clues to their exact nature. We first look at nonleptonic B decays 
involving the p-wave D s resonant states which we will denote by D s q, corresponding to the p-wave, 3q = \, 
+ state, and D* sl corresponding to the p-wave, j q = \, 1 + state. In B factories that do not produce the 
B s mesons the D s p-wave states cannot be directly produced via the weak current involving the b quark 
but they can only be produced through the sc current in the weak decay effective Hamiltonian. It was 
suggested in Ref. [16, 17] that these theoretically expected broad states may be discovered through the 
three body decays B — > D^D^K decays, where D^> refer to D or D*, if these states are above the 
D^K threshold. These three body decays can also be used to measure both sin 2(3 and cos 2(3 [18, 19, 16]. 
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In hadron B factories the D s resonant states can be produced directly from the weak decay of the b quark 
in the B s meson. 

In this work we concentrate on non leptonic decays of the type B — > £)WD s (2317) and B — > Z)Ml) s (2460), 
which are accessible at current B factories, and we also study nonleptonic decays of the types B s — > 
D S (2317)M and B s — ► D S (2460)M, where M is the meson formed by the emitted W. These latter decays 
can be studied at hadron B factories. Our purpose here is to explore what additional information about 
the structure and the properties of the new D s states can be obtained from these nonleptonic decays. 

2 Nonleptonic Decay 

Let us first assume that we can identify the the newly discovered states D s (2317) with D s q and D s (2460) 
with D* sl . In the Standard Model (SM) the amplitudes for B — > D^D s0 (D* 1 ), are generated by the 
following effective Hamiltonian [20]: 

(-i 10 

Hl f} = ^[^VftfaO^ + caO^ (1) 

v 1 i=3 
where the superscript u, c, t indicates the internal quark, / can be u or c quark, q can be either a d or a 
s quark depending on whether the decay is a AS = or AS = — 1 process. The operators Of are defined 
as [21] 

0\ } = WpLfpfp-fLba , O q 2f = q lfl Lff^Lb, 

Of i5 = q^Lbq'^L(R)q' , Of 6 = q allM Lb^ pl ^L{R)q' a , (2) 

3 3 
°7,9 = ■^.ln L be q >rfl IJ 'R{L)q' , O| 10 = -q^^LbpCq^^^R^q^ , 

where -R(£) = 1 ±75, and g' is summed over all flavors except t. Oi/,2/ are the current-current operators 
that represent tree level processes. O3-6 are the strong gluon induced penguin operators, and operators 
C>7_io are due to 7 and Z exchange (electroweak penguins), and "box" diagrams at loop level. The values 
of the Wilson coefficients can be found in Ref. [20]. 
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In the factorization assumption the amplitude for B — ► D s q(D* 1 ), can now be written as 

M = M 1 +M 2 (3) 



where 



where 



Mi = < D s0 (D* sl )\ s - 7 ^(l- 7 5 ) c | 0>< j dW|c 7 m (1-7 5 )M^> 

M 2 = ^X 2 <D s0 (D* sl )\s(l+ 1 5 )c\0><D^\c(l- 1 5 )b\B> (4) 

(5) 



^2 = -2 (-^-B 5 + £ 6 + + B 8 ) (6) 



We have defined 



with 



Bi = ~ £ c ^ (7) 



n = (8) 

In the above equations iV c represents the number of colors. To simplify matters we neglect the small 
penguin contributions and so as a first approximation we will neglect Mi- The currents involving the heavy 
b and c quarks, =< D|c7^(l — ^)b\B(p) > and =< D*(ei)\c^(l — j 5 )b\B(p) > can be expressed 
in terms of form factors [22]. In the heavy quark limit the various form factors are related to a universal 
Isgur-Wise function £(v-vi) where v and v\ are the four velocities of the B and the mesons. One can 
therefore write, 

J£ = V™bV™d£( v ■ v i) K + <] (9) 

and 

j£. = ^^i{v-v 1 )\-ie^e\ v v a v w + v^e\-v-e^{v-v l + l)\ (10) 



The matrix elements < D s q \ S7 M (1 — 7 5 ) c | > and < D*± \ S7 M (1 — 7 5 ) c | > are written in terms of 
the decay constants that are defined as 

< D s0 (P)\ s7m(1-7 5 )c|0> = if Da0 P^ 
<D* sl (P,e 2 )\s^(l- 1 5 )c\0> = M D *J Dti s*^ (11) 

We can now define the following ratios 

BR[B D£ s0 ] 



Rdo 
Rd*o = 
Rdi 



BR[B DD S ] 
£?i?[ff -> £)*D 80 ] 

" BR[B -> DDf\ 



BR[B -> D*D* al ] 

Rd * 1 ~ BR[B D*D*] (12) 
Let us focus on the ratio Rdo which within factorization and the heavy quark limit can be written as 

Rdo = l^l 2 (13) 

JD S 

where we have neglected phase space ( and other) effects that are subleading in the heavy quark expansion. 
Similarly we have 

Rdi = |^| 2 (14) 

JD* S 

Now in the heavy quark limit fo s0 = /d* s1 an d Jd s = Jd* and so one would predict Rdo ~ Rdi- 
There have been various estimates of the decay constant /d s0 in quark models [23] and in QCD sum rule 
calculations [24]; these typically find the p-wave , j q = \ states to have the similar decay constants as the 
ground state mesons. We therefore expect fo s0 ~ /d s giving in addition to the heavy quark predictions 

Rdo ~ Rdi ~ 1 (15) 

Experimentally Belle measures [11] 

BR[B -► DD S (2317)]BR[D S (2317) -► D s tt°] = (9.9±|;| ± 3.0) x 10" 4 (16) 



The dominant decay of the D s (2317) is expected to be through the D s ir mode [25, 26] and so 

BR[D -> DD 3 (2317)] « 1(T 3 (17) 
Now using the measured branching ratio [12] 

BR[B + -»• D°Df] = (1.3 ± 0.4) x W' 2 

BR[B d -► £>"£>+] = (8 ± 3) x 10~ 3 (18) 
one obtains a combined branching ratio 

BR[B — > DD S ] « 10~ 2 (19) 

This leads to i?oo ~ ^ ( Or >/A,o ~ s/dJ which is a factor 10 smaller then theoretical expectations. 
There are a few possible explanations that can be put forward to explain this discrepancy between exper- 
iment and theoretical expectation and we will consider them now. 

It is possible that the estimate of the decay constants of the p-wave, j q = \ states in the various models 
are incorrect just like the mass predictions of these states are incorrect. This would require a major revision 
of model calculations that predict the properties of these states. From the experimental data we have seen 
that f Ds0 ~ \f Ds which gives, using f Ds0 = f D ^, 

Rdi ~ ^ (20) 

To check this we note that experimentally Belle measures [11] 

BR[B -> DD S (2A60)]BR[D S (2A60)} L> s tt ] = (25.8±J;g ± 7.7) x 10~ 4 
BR[B -► DD S (2A60)]BR[D S (2460)] -► L> s7 ] = (5.3±i;| ± 1.6) x 10" 4 (21) 

Taking the central values we find 

BR[B + -> T?D a (246Q)] < 31.1 x 10~ 4 (22) 
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Using the measured branching ratio [12] 



BR[B+ -► WD+*] 



(9 ± 4) x 10~ 3 



BR[B d - D-D+*] 



(1.0 ±0.5) x 10 



-2 



(23) 



one can obtain, using the measured central values 



BR[B ->■ DL>*] 



10 



,-2 



(24) 



This then leads to Rdi 



| which is in disagreement with Eq. 15 and Eq. 20. 



One might argue that factorization is not applicable to B — > £)(*)£)(*) decays. However recent analysis 



in Ref. [27] find that factorization works well for these decays. Moreover the quantities in Eq. 12 are 
ratios of nonleptonic decay amplitudes and so nonfactorizable effects may cancel. So what one really 
requires is significantly different nonfactorizable corrections between decays with the p-wave states in the 
final state and decays with the ground state mesons in the final state. It is possible that the discrepancies 
between experiments and theory may arise from a combination of incorrect model prediction of p-wave 
state properties and nonfactorizable effects. 

3 Nonleptonic decays involving B s decays 

Another test of the the nature of the newly discovered D s states that does not rely on factorization or 
heavy quark symmetry involves the B s mesons. As we indicated earlier, with the B s meson, the p-wave 
D s states can be produced via the weak current involving the b quark. We can now consider decays of the 
type B s -► L> S (2317)(L> S (2460))M where M = ir,p,K etc. With the identification of D s (2317) (D s (2460)) 
as the p-wave states these decays are the same as B s — ► D s q(D* 1 )M. 

One can now construct the ratios 




D s (2317) 



BR[B S -► L> S (2317)M] 
BR[B d D d0 M] 

BR[B S -► L> S (2460)M] 
BR[B d -+ D* dl M] 
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= BR[B S -> D S M] 
Ds BR[B d D^M] 
_ -> D* S M] 

Dl BR[B d D*M] 1 J 

Now in the S77(3) limit all the ratios are unity. Moreover the ratio of ratios r = D ^ 317) and n = p ^f 0) 
are expected to have smaller flavour symmetry violations and hence smaller deviations from unity, as SU (3) 
breaking effects in the ratios may cancel [28]. Hence any large deviation of T Ds ^ 2317 ^ and r Ds ( 2 46o) from 
unity would be inconsistent with the j q = \ p-wave interpretation of the new D s states. Note that the 
further assumption of factorization leads to T D j 2 zi7) ~ ^£> 3 (2460) an d ^d s ~ Tp* in the heavy quark limit. 

As indicated earlier, among various other suggestions for the nature of the new D s states is the idea 
that these states may be D^K molecules. There are no serious models of such meson molecules that one 
can use to calculate nonleptonic decays involving these states. Here we will attempt a rough qualitative 
estimate of nonleptonic decay rates assuming that the _D S (2317) and A, (2460) states are really a DK 
molecule and a D*K molecule respectively. Consider the nonleptonic decay B — > DD S (2317). We assume 
that the decay proceeds through two stages: the first stage is the decay B — > DDK, followed by the state 
D(p 2 )K(pk) forming the molecule D s (2317) with the probability given by /(p2,Pk) so that 

dY{B - £>£> a (2317)) = 1 ^ \A(B -> D(p 1 )D(p 2 )K(p K ))\ 2 f(p 2 ,p K )dE K dE 2 (26) 

Without a model for f(p2,Px) we cannot make predictions but nonetheless it is useful to define the average 
probability function / as 

? I \MB -> D( Pl )D(p 2 )K(p K ))\ 2 f(p 2 ,p K )dE K dE 2 

1 J\A(B - D( Pl )D(p 2 )K(p K ))\2dE K dE 2 ( > 

Hence we have 

BR(B+ D i °* ) D s (2317) + ) = BR(B + D {0 * ] D + K ) x J 

BR(B° -► D^D S (2317) + ) = BR(B° -► D^D°K + ) x J (28) 
We can define a similar function /* and the average /* for nonleptonic decays involving the -D s (2460) and 



so 



BR{B + -► £> (0 * ) .D i (2460) + ) = BR{B + D {G * ] D + *K Q ) x /* 
Si?(B -► £>(-*) D s (2460)+) = BR(B° ^ D ( --*^D°*K + ) xj* 



We can consider the ratios 



z+ = ffit!(i? + ^ 7^(2460)+) z+ „ = 5fl(g+^^°*£) s (2460)+) 
re " -^7J°D S (2317)+)' reS BR{B+ -^^* D s (2317)+) 

z o = BR(B ^ D-D 8 {2460)+) ^ = BR(B° -> D-*D S (2A60)+) 
res BR(B° ^D~D s (2317)+y res BR(B° -> D~* D s (2317)+) 

z+ BR{B+ -> D°D + *K°) z+ * BR(B + -> D°*D+*K°) 

o _ Bi?(5° -► D-D°*K + ) ^ __ BR(B° -► D~*D°*K + ) 
3 - b ° dv ~ BR(B°^D-D°K+) 3 - b ° dy ~ BR(B°^D-*D°K+) 



which are related as 



7 7 

7+ 7+ £ y+* 7+* £_ 

^res — ^3-body J ' ^res ~ ^3-body J 



yO yO f yO* yO* f 

^res — A 3-body^ji A res ~ ^3-body^j 



Using the measured three body branching ratios [29] 
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BR(B + -► D D + K ) = (0.18 ± 0.07 ± 0.04) x 10 

BR(B° — > D~D°K + ) = (0.17 ± 0.03 ± 0.03) x 10~ 2 

BR(B+ D°*D + K°) = (0.41±g;i| ± 0.08) x 10~ 2 

5i?(£° -► D-*D°K+) = (0.3llg;^ ± 0.04) x 10~ 2 

Si?(5 + -► lfD + *K ) = (0.52±g;Jg ± 0.07) x 10~ 2 

BR(B° — > D~D*°K + ) = (0.46 ± 0.07 ± 0.07) x 10" 2 

-► D°*D + *K°) = (0.78±g;|f ± 0.14) x 10~ 2 

BR(B° — ► D~* D°* K + ) = (1.18 ± 0.10 ± 0.17) x 10~ 2 
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which are proportional to either 1 — / or 1 — /' and assuming /* « /* allows one to obtain, with the 
central values of the measurements, 

Ztbody = 2 - 89 (33) 
which can be compared to Z+ s = 3.14 from Eq. 17 and Eq. 22. If fact the prediction Z+ es ~ 3 , Z^* s ~ 3, 
%res ~ 3 and Z®* s ~ 3 are consistent within the errors for the three body branching ratios in Z^^^y, 
Z t-body Z 3-bod y and Z 3-body We also obtain / ~ /* ~ 0.3 , from Eq. 17 and Eq. 32 which indicates that 
a sizeable fraction of the D^K state form molecules. 

Finally we can extend this model also to the case where the D s resonance is produced via the weak 
current containing the b quark in B s decays. Consider the decays B s — ► D S (2317)M where M is the 
emitted meson. The form factor for B — ► D s (2317) transition can then be related to B — > DK transition. 
In other words, we can write 

BR[B S D s (2317)Af] = BR[B S ^ DKM]f' (34) 

where 

r J |A(g ^ Zj(p 2 )^(p y )M( Pl ))| 2 /(p 2 ,p x )^^2 , , 

7 j\A{B ^ D{p 2 )K{p K )M{ Pl ))\HE K dE 2 { ' 

— 

We can similarly define / as 

Bi?[5 s -► L> S (2460)M] = 5i2[fi s -► D*KM]f*' (36) 

Note that the ratios T Ds ( 2 3i7) an d ^d s (2460) (Eq- 25) in the molecular model are no longer equal to unity 
in the SU(3) limit since that depended on the identification of these states as p-wave states. Therefore 
the measurement of these ratios can provide useful information on the nature of the D s (2317) and the 
A, (2460) states. 

4 Summary and Conclusions 

In summary, in this work, we have considered the nonleptonic B decays B -> D S (2317)(D S (2460)), 
involving the newly discovered D s (2317) and the D s (2460) states. We have discussed the implication of 
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the measured nonleptonic decays for the properties and the nature of these states. If these states are the 
p-wave multiplet with the light degrees of freedom in the j q = \ state, then we find that experiments 
indicate disagreement with model calculation of their properties and/or breakdown of the factorization 
assumption. We have suggested further tests involving nonleptonic B s meson decays, that do not assume 
factorization but assumes SU(3) flavour symmetry, that can further shed light on the true nature of these 
newly discovered states. Finally, we have also proposed a model to calculate the two body nonleptonic 
decays B -► £>(*) D S (2317)(L> S (2460)), assuming that the the D s (2317) and L> s (2460) are DK and D*K 
molecules. The model relates these two body nonleptonic decays to the three body B decays of the type 
B -> DWdWR. 
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